Dissipation-induced quantum phase transition in a quantum box 
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In a recent work, Le Hur has shown that dissipative couphng to gate electrodes may play an 
important role in a quantum box near its degeneracy point [K. Le Hur, Phys. Rev. Lett. 92, 
196804 (2004)]: While quantum fluctuations of the charge of the dot tend to round Coulomb blockade 
charging steps of the box, strong enough dissipation suppresses these fluctuations and leads to the 
reappearance of sharp charging steps. In the present paper we study this quantum phase transition 
in detail using bosonization and numerical renormalization group methods in the limit of vanishing 
level spacing. 



PACS numbers: 75.20.Hr, 71.27.+a, 72.15.Qm 

I. INTRODUCTION 

Coulomb blockade is one of the most studied and 
most basic correlation effects that occurs in mesoscopic 
devicesi This phenomenon appears in small structures 
(quantum dots or metallic islands) weakly connected to 
the rest of the world. The capacitance of these small de- 
vices can be very small and therefore the charging energy 
Ec, corresponding to the energy cost of putting an extra 
electron on the device, can become large. If this charging 
energy is larger than the measurement temperature then 
correlation effects can become important, thermal fluc- 
tuations of the charge of the island become small, and 
transport through the device is typically suppressed as 
we lower the temperature. 

In addition to the above-mentioned thermal fluctua- 
tions. Coulomb blockade can also be lifted by quantum 
fluctuations. There is two important known mechanisms 
that lead to quantum fluctuations: (a) Depending on the 
applied gate voltages, two charging states of the quantum 
dot can be degenerate, and electrons can thus gain kinetic 
energy by hybridizing with states in the leads, accompa- 
nied by strong quantum fluctuations of the charge?2i^ (b) 
Quantum fluctuations can play an important role at tem- 
peratures below the single particle level spacing of the dot 
too, if the dot contains an odd number of electrons or if 
the ground state of the isolated dot is degenerate. In this 
case the system gets rid of the residual entropy associated 
with the spin degeneracy of the isolated dot through the 
Kondo effect (4 which consists of binding conduction elec- 
trons in the attached metallic leads antiferromagnetically 
to the spin of the dot.^'*^ 

The simplest device where charge fluctuations play a 
dominant role is the single electron box (SEB) schemat- 
ically shown in Fig. ^ In this device one couples an iso- 
lated quantum box capacitively to a gate electrode and to 
an external lead through a tunnel junction. The isolated 
box can be very well described by the following simple 
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FIG. 1: Sketch of the single electron box. Fig. b. shows 
the top view of the regime where electrons move in the the 
two-dimensional electron gas. Black areas indicate various 
electrodes necessary to shape the electron gas. 
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where dj^ ^ creates an electron on the dot with spin a 
in a single particle state of energy e„. The second term 
accounts for the electron-electron interaction with C the 
total capacitance of the box, and N = ^ : (il^^dn,^ ■ 
the number of extra electrons on the island. In Eq. 
Ug = VgCg/e stauds the dimensionless gate voltage ex- 
pressed in terms of the gate voltage Vg , the gate capaci- 
tance, Cg, and the electron charge e. Clearly, the number 
of electrons on an isolated box jumps by one at approxi- 
mately half-integer values of Ug as we increase Ug^ 

In reality, however, the box is not fully isolated from 
the rest of the world, and in the vicinity of these degener- 
acy points quantum fluctuations to the leads need to be 
considered. In the following we shall restrict our consid- 
erations to the particularly interesting regime where the 
temperature is much larger than the typical single parti- 
cle level spacing A on the box, but is also much smaller 
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FIG. 2: The number of electrons on a single electron box as 
a function of the dimensionless gate voltage, (a) The sudden 
jumps of an isolated box become smeared out due to quantum 
fluctuation as soon as we couple the box to a lead, (b) Large 
enough dissipation restores the discrete jumps. 



than the charging energy Ec — e^/2C, 
A < T < £;c . 



(2) 



Furthermore, we shaU focus our attention to the vicin- 
ity of the above degeneracy points, Ug w half-integer. 
As shown by Matveevp in this regime quantum fluctua- 
tions of the charge on the dot can be treated by mapping 
the Hamiltonian to that of the anisotropic two-channel 
Kondo problem. The finite jumps in (TV) are rounded 
by these quantum fluctuations, and all what is left from 
them in the A ^ limit is a logarithmic divergence in 
the slope of {N){ng) at the degeneracy points rig « half- 
integer (see Fig. |21a). 

However, as shown recently by Le Hur, this picture is 
not complete:^ In Eq. one typically assumes that the 
gate voltage is constant. In reality, however, this gate 
voltage represents electromagnetic degrees of freedom. 
The fluctuations of these electromagnetic degrees of free- 
dom must be considered and typically result in Ohmic 
dissipation, characterized by a dimensionless dissipation 
strength a.i^^ The specific value of this dissipation de- 
pends on the resistance Rg of the gate electrode and the 
gate capacitance, and is proportional to 
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(3) 



with Rq — the quantum resistance unit. As shown 
by Le Hur, this additional dissipative coupling leads to 
a quantum phase transition: Below a critical dissipation 
strength, ac, dissipation is ultimately irrelevant, and the 
scenario of Matveev holds with somewhat modified pa- 
rameters. If, however, a is large enough then dissipa- 
tion dominates, charge fluctuations are suppressed, and 
{N)(ng) displays finite jumps at half- integer values of rig 
at T = temperature, as shown in Fig. [21b. ^ This crit- 
ical value of dissipation depends on the dimensionless 
conductance g of the tunnel junction, ac — oic(g)- 

We have to remark here that the above considera- 
tions hold only in the limit of vanishing level spacing, 
A — > 0. In reality, level spacing on the dot is finite, and 
for semiconducting devices the ratio Ec / ^ cannot be 



much larger than about lOO.-^-^- Thus the level spacing 
A ultimately provides a finite cut-off similar to a finite 
temperature, and the slope of {N){ng) is expected to re- 
main finite even at T = temperature. Although it is not 
clear what happens below the level spacing, in a real sys- 
tem the phase transition above is thus probably reduced 
to a cross-over (see our conclusions too). 

In the present paper we shall investigate this phase 
transition employing non-perturbative methods such as 
Abelian bosonizationi^ and numerical renormalization 
group to a somewhat generalized version of the quan- 
tum box modeli22i2i These methods enable us to obtain 
the full phase diagram of the quantum box. 

Using the bosonization approach we are able to con- 
struct the non-perturbative scaling equations that de- 
scribe the above quantum phase transition. At the de- 
generacy point, we find a Kosterlitz-Thouless quantum 
phase transition, in agreement with the weak coupling 
result of Le HuriS We also determine the scaling of the 
Kondo temperature Tk associated with the charge fluc- 
tuations non-perturbatively. Apart from an overall pref- 
actor, this scale vanishes as 



TA'(a) ~ Ec exp|- 



}, (4) 



where g — j^G is the dimensionless conductance of the 
junction. Thus Tk approaches zero very fast but con- 
tinuously as one approaches the critical value ac from 
below. 

The phase transition also manifests itself in the T ~ Q 
temperature linear conductance of the tunnel junction, 
which displays a first order phase transition: 



g{T = 0,a) = ^^("c - a) . 



(5) 



We also study the above phase transition using the 
powerful method of numerical renormalization group 
(NRG). Since the calculation for fermions with spin is 
computationally extremely demanding, we will have to 
restrict ourself to the case of spinless fermions. The phase 
transition within this method simply appears as a change 
in the finite size spectrum. The NRG method also en- 
ables us to determine the shape of the step in {N){ng) 
and Tk non-perturbatively, and compute various spectral 
functions. Our results agree very well with the scenario 
of Ref. 9 and the bosonization results. 

The paper is structured as follows: In Section ITU we 
introduce our basic notations and show how to map the 
problem to a Kondo model. Section IIIII is devoted to 
the construction of the non-perturbative scaling equa- 
tions and their analysis. The numerical renormalization 
group results are presented in Section IIVI Finally, we 
summarize our results and discuss future perspectives in 
Section El 
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II. HAMILTONIAN 

We shall assume that couphng between the lead and 
the dot is small, so that charge fluctuations can be de- 
scribed within the tunneling approximation 



= t 



E 



(4,. 



+ h.c.) , 



(6) 



where the operators cl ^ create conduction electrons in 
the lead with energy e and spin a, 



(7) 



In Eq. © we assumed that there is only a single tunneling 
mode that is coupled to the dot electrons so that the 
tunneling matrix elements can be replaced by a simple 
number. The situation is much more complex in case of 
several tunneling modes and shall not be discussed here 
(see Ref.0)- 

Since all half-integer values of ng are essentially equiv- 
alent, to obtain Matveev's mapping, we shall focus to 
the vicinity of the point Ug ~ 1/2. The important quan- 
tum fluctuations in this regime are dominated by the two 
charging states |iV = 0) and \N = 1) of the dot. Other 
charging states of the dot are separated at least by an 
energy ~ Eq, and in the small temperature/tunneling 
regime their effect is only to slightly renormalize the val- 
ues of t and Ec- We can thus restrict our considerations 
to the charging states = and = 1 , and keep track 
of charge fluctuations of the box by introducing an or- 
bital pseudospin T, the states = (1 - N)/2 = ±1/2 
corresponding to states with = 0, 1. Note that we still 
have to keep track of internal electron-hole excitations of 
the box, since our energy scale is larger than the level 
spacing. 

To rewrite the tunneling Hamiltonian in a more sug- 
gestive way we introduce the new fields normalized by 
the density of states in the box and in the lead, gho^ and 
ffload, respectively 



1 



1 



and organize them into a four component spinor 



(9) 



The normalization in Eq. ||SJ| has been chosen is such a 
way that the imaginary time propagator of the field 

satisfies (Tt-V't(t(t)V'I'ct' (*-*)) = 5a^a'^T,T'/T ■ In this lan- 
guage we can rewrite the tunneling part of Hamiltonian 
as 



(T+V4T-^. + h.c.) 



(10) 



where the operator just fiips the orbital spin r of the 
field ^l^ra, and the j'j^ is a dimensionless coupling propor- 
tional to the tunneling, = 2iY/pb^J^^^ijr^. Thus j\ is 



directly related to the dimensionless conductance g of the 
tunnel junction 



G 7r2 , 



(11) 



with Gq = 2e^//i, the quantum conductance unit. 
Clearly, Eq. (|l()|l is just the Hamiltonian of an anisotropic 
two-channel Kondo modelfSiii the orbital spins T and r 
playing the role of the spins of the original two-channel 
Kondo model, and the electron spin a providing the silent 
channel index. The presence of this additional channel 
index makes the physics of the two-channel Kondo model 
entirely different from that of the single channel Kondo 
problem, and leads to non-Fermi liquid properties, such 
as the logarithmically divergent capacitance in our case. 
We emphasize again, that the above mapping holds only 
in the regime A < T, w, ... < Eq, where the level spacing 
of the box can be neglected. 

Deviations from the point Ug — 1/2 appear in this 
language as a local field acting on the orbital spin of the 
box 



2C(^ 



-BT^ 



(12) 



where the effective magnetic field is simply given by 
B = Ec{^ — 2ng). It is clear from this expression that 
fiuctuations of the gate voltage, 5Vg will result in a dis- 
sipative coupling of the following form 



diss 



(13) 



where for an Ohmic heat bath the bosonic field ip decays 



(T,^(r)^(0)) = ^ 



(14) 



The dimensionless coupling A above is related to the 
usual couphng a as a = A^/2. The specific value of A 
depends on the properties of gate electrodes and is ap- 
proximately given by Eq. Q . Thus for very resistive gate 
electrodes, a can be large, while if the gate electrodes are 
good conductors, then a is small. 

Eqs. fTn|l . (fT^ . and (fnji constitute the effective Hamil- 
tonian that describes the physics of the dot in the vicinity 
of the degeneracy point in the presence of dissipation. In 
the rest of the paper we shall, however, include an addi- 
tional term in the Hamiltonian of the form 



(15) 



too. This term describes charging state dependent scat- 
tering of the conduction electrons. While this term is 
typically small, it is generated under the renormalization 
group procedure. Furthermore, in other physical systems 
such as a Kondo spin coupled to fluctuations of quan- 
tum critical spin degrees of freedom, e.g.^ this term is 
present in the bare Hamiltonian as well. 
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III. BOSONIZATION 

The non-perturbative method of AbeUan bosonization 
can be very efficiently used to construct the phase dia- 
gram of our model, and allows us to treat both the dissi- 
pative coupling A and the coupling jz non-perturbatively. 
For the sake of simplicity, we shall first discuss the some- 
what simpler case of spinless Fermions and then restore 
spin indices. Throughout this Section we shall use the 
bosonization scheme of Ref. UM 



A. Bosonization for spinless Fermions 

For spinless fermions the field ^ has only two orbital 
indices, r = ±. Correspondingly, we have to introduce 
two left-moving bosonic fields $± to represent ■(/'±r^ 



\/a 



(16) 



where F± denotes the Klein factor, and a is a short dis- 
tance cut-off ~ vf/Ec- In terms of these bosonic fields 
the interaction part of the Hamiltonian reads 



n/2 



1 



(17) 

+ h.c.), (18) 



where for the sake of simplicity we suppressed Klein 
factors and introduced the orbital field, $t — (^+ — 
^-)/V2. Note that the charge field $c = {<^+ + <^-)/V2 
completely decouples from the orbital spin of the box. 
Particular care is needed while deriving Eq. H17|l : Intro- 
ducing the bosonization cut-off scheme renormalizes the 
coupling constants, jz ^ jz- 



'iSz 4 fnjz 
jz = = -atan — 

TT TT V 4 



(19) 



with 5z the phase shift generated by jz- Representing 
the field </? as the derivative of another bosonic field <E> we 
can finally rewrite the dissipative term as 



(20) 



In this bosonizcd form the non-interacting part of the 
Hamiltonian becomes 



•J i-r r~< 



(21) 



where the first two terms represent electron-hole excita- 
tions on the dot and in the lead, and the last term gives 
the energy of excitations responsible for dissipation. 

To obtain the phase diagram of our model we carried 
out a non-perturbative scaling analysis. As a first step we 



eliminated the terms proportional to A and jz by applying 
a unitary transformation 

V = exp{i Tz (A$(0) + i^$7.(0)) } . (22) 

This unitary transformation brings iJpcrp to the following 
form 



+ ^A<E.(0)+^(^-^/2)<i>T(0) 



+ h.c. 



(23) 



The scaling dimension of j^ can be read out of this 
equation and one obtains 



dl 



j± + 0(j1) 



(24) 



where a = A^/2 is the dissipation strength, and / = In a 
is the scaling variable. The first term in this equation is 
just the famous non-perturbative equation of Yuval and 
Andersonpi^ The second term in this scaling equation just 
means that dissipative coupling to charge fiuctuations on 
the gate tends to suppress quantum fluctuations of the 
charge of the dot. It is this term that drives the phase 
transition. 

The scaling equations of a and Sz can be simply ob- 
tained by taking the operator product expansion of the 
term j± with itself. In this way we obtain: 



/4^\ ^ 'Mz 



dl\TT 

da 2 



2(5, 



(25) 
(26) 



where terms of order j'j_ have been dropped. The 
first of these equations has been derived by Yuval and 
Anderson, and shows that j± tends to boost up Sz and 
generate a Kondo effect. The second equation shows that 
quantum fluctuations suppress dissipation effects. There 
is thus a competition between j±_ and a, which mutually 
tend to suppress each-other. 

Depending on the model parameters, we have to dis- 
tinguish two regimes: In the small dissipation regime dis- 
sipation is irrelevant, a scales to zero, and charge fluctu- 
ations generated by j± ~ t start to dominate at an en- 
ergy scale Tk, the so-called Kondo energy. This Kondo 
fixed point corresponds to 6z = tt/2 and j± 0{l)m^ 
It is easy to show that the dissipative coupling is indeed 
irrelevant at this Kondo fixed point: This follows from 
the observation that the Kondo fixed point is described 
by a Fermi liquid theory»4 As a consequence, the orbital 
susceptibility is finite at zero temperature, xt ~ ^/Tk, 
and correspondingly the imaginary time correlation func- 
tion of Tz decays as (T^r^(T)r^(O)) ~ l/r^. Therefore 
the correlation function {TrTz{T)ip{T)Tz{0)ip{0)) decays 
as 1/t*, implying that the scahng dimension of the dis- 
sipation is —1 at the Kondo fixed point, and is therefore 
irrelevant lii The capacitance of the dot for these small 
values of a is simply proportional to the orbital suscep- 
tibility, C^XT- I/Tk- 
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If, on the other hand, the bare value of a is larger 
than a critical value, acO'_L,(52), then quantum fluctua- 
tions are suppressed, and jj_ scales to implying that 
charge fluctuations are suppressed. This phase can be 
called localized in the sense that if we prepare the box in 
the state iV = at the degeneracy point, then the ex- 
pectation value of N will never approach its equilibrium 
value 1/2, 



while the other two scaling equations read 



lini (N(t)) < 1/2 

t — >C30 



This behavior is analogous to that of the dissipative two- 
state system.^" 

It is easy to determine this critical value of a in the 
j±_ limit from our scaling equations: 



"cOx ^ 0) = 4- 



1 



(27) 



To obtain the phase diagram in the vicinity of this line, 
we rewrite the scaling equations as 



dl 
dr 

'dl ~ 



j1 



(28) 
(29) 

(30) 



where we introduced the variable 

TT V TT / 

measuring the distance from the line Eq. H27|) . These 
equations clearly describe a Kosterlitz-Thouless phase 
transitionii^ Solving these equations we find that for 
small values of j± 



bl|+4^(l-^) 



(31) 



and the corresponding Kondo scale vanishes as one ap- 
proaches a as 



Tk ~ Ec exp|- 



(32) 



It is important to remark at this point that the above 
scaling equations describe the nature of the phase tran- 
sition even for large values of j± . The reason is that in 
the vicinity of the quantum phase transition all scaling 
trajectories approach the critical line defined by j± = 
and a = 4^[1 - Thus the initial part of the scahng 
just renormalizes the effective value of the various cou- 
plings that enter Eq. (|40|l . but the overall picture does 
not change. 



B. Fermions with spin 

The considerations of the previous subsection can be 
easily generalized to the case of Fermions with spin. The 
scaling equation of j± then becomes 



dl 



. TT 



a 



(33) 



d_/45, 

dl\ TT 

da 



= i- 



4(5, 



3l 



(34) 
(35) 



In the absence of dissipation, these scaling equa- 
tions have been derived first by Vladar, Zimanyi and 
Zawadowski.^^ Similar to the spinless case, two regimes 
can be distinguished!^ For small dissipations a is irrel- 
evant, and at the degeneracy point of the dot the cou- 
plings j± and 5z flow to the two-channel Kondo fixed 
point, Sz = 7r/4, j± 0{1). The physics is, however, 
not that of a Fermi liquid due to the presence of electron 
spin: the orbital susceptibility {i.e. the capacitance) of 
the dot diverges logarithmically below the Kondo tem- 
perature, XT ^ \n{TK/T)/TK^ This logarithmic di- 
vergence indicates that the deviation from the degener- 
acy point rig = 1/2 represent a relevant perturbation of 
scaling dimension 1/2, and as a consequence, the slope of 
(TV) (rig) diverges logarithmically at T = temperature 
as one approaches = 1 /2 (if we neglect the finite level 
spacing on the dot). 

It is easy to check that dissipation is irrelevant at this 
two-channel Kondo fixed point too: The logarithmic di- 
vergence of the orbital susceptibility corresponds to an 
asymptotic decay {TrTz{T)Tz{0)) ^ 1/r. Consequently, 
the correlation function {T rTz{T)ip{T)Tz{0)(p{0)) decays 
as 1/t^ at the two-channel Kondo fixed point, implying 
that the scaling dimension of the dissipation is —1/2 at 
the Kondo fixed point, and is therefore irrelevant lii 

If, however, a is large enough then again, the dissipa- 
tion drives j± to zero, and we recover the localized phase 
with jumps in the function {N)(ng). Similar to the spin- 
less case, we can determine the boundary between these 
two phases by introducing the variable 



TT V TT J 



(36) 



(37) 



and rewriting the scaling equations in terms of it. Keep- 
ing only the lowest order terms in r and j± we obtain 
again a set of two coupled differential equations describ- 
ing a Kosterlitz-Thouless phase transition: 

dj±_ 
dl 

For small values of j± the phase boundary is thus given 
by 

af"«jl + 4^(l-— ) , (39) 

TT \ TT / 

and the Kondo scale vanishes as one approaches etc from 
below as 



Tk w Ec cxp\ - 



{- 



y'2|j^|(l-ar")i/Var"-a 



} 



(40) 
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IV. NUMERICAL RENORMALIZATION 
GROUP CALCULATION 

In order to obtain results for different physical quan- 
tities we performed Wilson's numerical renormaliza- 
tion group method (NRG) on the modeli22i2ii22i2^ This 
method is especially suitable for quantum impurity prob- 
lems and can be used to determine dynamical as well 
as thermodynamic properties with high accuracy. How- 
ever, performing such a calculation on the Bose-Fermi 
Kondo model, a conceptual problem arises: Originally, 
NRG was designed to treat fermionic problems, where 
the method is based on a logarithmic discretization of 
the conduction band and then a mapping of it onto a 
semi-infinite chain. In the present model we need to treat 
both fermionic and bosonic baths at the same time. Even 
though the NRG method has recently been extended to 
bosonic systems^^^ the bosonic calculation needs con- 
siderable computational effort. We therefore decided to 
treat this problem in a different way. 

Fortunately, in the present model the bosonic field (p 
is Ohmic. This implies that we can represent this field 
by coupling Tz to the total density of F fermionic fields 
in the following way: 

F 

i?di.s = |T^5]*I*.. (41) 

1=1 

Here the fields denote some fermion fields normalized 
as (TT-5'i(T)^'j(0)) — Sij/r. To make this mapping com- 
plete, we have to find the relation between q and the 
dissipation strength a. It is a relatively simple matter to 
show that these are connected as 

a = 2F(^^^\ (42) 

with 5g ~ atan(7r(7/4) the phase shift associated with the 
potential scattering g. 

The advantage of this trick is that we can use just 
fermionic fields in our calculations. On the other hand, 
this trick also implies that all along the calculations, we 
carry along F—1 bosonic fields which are decoupled from 
Tz- For large values of F we thus keep track of many 
redundant degrees of freedom. Furthermore, we cannot 
reach any value of a in this way. In our calculations, 
e.g., we used F = 2 fermion fields allowing to map out 
the phase diagram for < a < 1. 

It is clear from the previous subsection that to repre- 
sent the physical system considered in the present paper, 
one needs at least three bands of spinful fermions: one of 
them, ^'o- would be used to represent dissipative term in 
the Hamiltonian while two others are just the electronic 
degrees of freedom, tpra- This imphes that one needs a 
three-channel NRG code to attack the problem in its full 
glory. Unfortunately, we do not have too many symme- 
tries simplifying our problem: Throughout the calcula- 
tion we used an SU(2) symmetry associated with particle- 
hole symmetry and a U(l) symmetry associated with the 
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FIG. 3; (color online) Relation between the phase shift 5g 
extracted from the finite size spectrum and the fermionic cou- 
phng g. The continuous line shows the fit with Eq. (I43II . 

z-component of the orbital spin. As a consequence, we 
were able to study only the case of spinless fermions with 
sufficient accuracy. Fortunately, as shown in the previ- 
ous section, most features of the spinful case already ap- 
pear in the spinless version of the model. Furthermore, 
the spinless case is also of relevance: this is realized in 
a single electron box by applying an external magnetic 
field that lifts the electron-spin degeneracy and drives the 
system away from the unstable two-channel Kondo fixed 
point. 

The very first step is to calibrate the method, i.e., to 
determine the correspondence between the fermionic cou- 
pling g and a. In principle, these are related by Eq. H42(l . 
However, the NRG discretization renormalizes the effec- 
tive value of g. To do that, we have performed calcula- 
tions for j_L = = in the presence of a term AqT^, 
i.e., for the spin-boson problem, and extracted the phase 
shifts from the finite size spectra. In that way we found 
that the phase shifts can be expressed with very high 
accuracy as 

5g = atan(/(A) g) , (43) 

where A is the parameter of the logarithmic discretization 
used in NRG and / (A) is a factor close to unity that must 
be determined numerically. In particular, for A = 2 used 
throughout this paper we find /(A = 2) — 1.03. The 
phase shift above can then be plugged into Eq. H42|l to 
obtain the corresponding value of a. This procedure is 
illustrated in Fig|31 A similar expression gives the value 
of 5z as a function of jz . 

The most basic result one can obtain by NRG is the 
finite size spectrum, i.e., the rescaled lowest lying energy 
levels plotted against the iteration number. We show two 
typical spectra in Fig0| In the top panel of Fig^ the 
spectrum crosses over to that of the usual Kondo fixed 
point^° with uniform level spacing at around iV ~ 60. 
In the bottom panel, however, no such crossover occurs, 
at least not before the 80th iteration indicating that the 
Kondo effect is suppressed by the strong coupling to the 
bosonic bath which tends to suppress charge fluctuations 
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FIG. 4: (color online) Finite size spectra obtained by NRG for 
the two distinct regimes. The upper part shows the spectrum 
for a = 0.653, jz = 0.2, j± = 0.7, while the lower part stands 
for the same js but for a — 0.667. It is transparent that 
while on the left figure the levels crosses over to the well- 
known Kondo spectrum, on the right panel no such crossover 
occurs. 



and thus localizes the orbital spin of the dot. This ro- 
bust difference in the finite size spectra allows us to ex- 
tract the critical value of the bosonic coupling ac and 
map out the phase boundary for all values of j± and Sz ■ 
Numerically we say that we are in the localized (dissi- 
pation dominated) phase if the crossover to the Kondo 
spectrum does not occur within TV = 80 iterations. The 
critical surface corresponding to this phase transition is 
shown in Fig.©. 

The case jz = is of particular interest, since this 
corresponds to the charge fluctuations at the degeneracy 
point of the quantum dot. In agreement with Eq. ac 
increases approximately linearly with j± ^ t y/g see 
Fig. . In Fig. |H1 we also show the phase boundary as 
a function of (J^/tt. The results compare very well with 
Eq. PHl . 

The dissipation dependence of the Kondo temperature 
can be also easily extracted from the finite size spec- 
trum. In Fig. Owe show l/\n{Ec /Tk) as a function of a. 
According to Ea.( I32|) this should vanish as ~ a/qic — a 
as one approaches the phase boundary, implying an ex- 
tremely fast change in Tk as one approaches ac- Though 
it is very difhcult numerically to get to the real scaling 
regime, the obtained results are quite consistent with the 
analytical expression Ea.(l32|l. 

To compute the full shape of the {N){ng) curve, one 
has to study the model for gate voltages Ug ^ 1/2, which 
corresponds to a field _B = (1 — 2ng)Ec applied to the 



FIG. 5: (color online) The critical surface in the parameter 
space separating the Kondo regime from the localized regime. 
The heavy continuous line indicates the phase transition line 
for the case relevant for the single electron box jz = 0. The 
heavy dashed line indicates the critical line, j± = 0. The 
arrows show the scaling trajectories for a — a^. All these 
critical scaling tranjectories end up at the critical line. 
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FIG. 6: Top: (color online) The physically relevant jz = cut 
of the critical surface. Bottom: Phase boundary as a function 
of Sz/iT for two values of j±. These phase boundaries compare 
very well with the analytical expression Eq. I|81|l . 
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FIG. 7: (color online) Kondo temperature as a function of 
dissipation, for j± = 0.4 and jz — 0.1 
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FIG. 8: (color online) The expectation value of A'^ — 1/2 for 
jz = 0.2, j± = 0.8 and different values of a. For moderate 
values of the bosonic coupling the step is smeared out by 
quantum fluctuations of the charge of the dot. The step gets 
sharper and sharper as a approaches ac, and For a > ac a, 
jump appears in (N). 



orbital spin = (1 — N)/2 within our approach. The 
result is shown in Fig|S| Clearly, for a < ac the Kondo 
correlations smear out the charging step of the grain. 

The Kondo temperature decreases as we increase the 
coupling of the spin to the bosonic environment, and re- 
sults in the sharpening of the step in (N). If we increase 
the bosonic coupling further, the Kondo effect disappears 
and {N) jumps at Ug = 1/2. It is difficult to tell from 
the numerics if the jump remains finite as a approaches 
ac from above. 

However, since the phase transition is of Kosterlitz- 
Thouless type, one expects that the 'order parameter', 
i.e. the size of the jump approaches a constant as a 
approaches ac from the dissipation-dominated (localized) 
sideS. 

The NRG method also enables us to compute local cor- 
relation functions, directly related to measurable quanti- 
ties. To determine the scattering properties of the con- 
duction electrons, we determined the conduction elec- 
trons' T-matrix. In the present case, the imaginary part 
of this quantity is related to the renormalized conduc- 
tance of the box-lead junction. The computation of this 
T-matrix can be reduced to the determination of the so- 
called 'composite Fermions' spectral function that can be 
computed by NRG similar to the single channel Kondo 
problemf2L2£ 

ImT^(u;) = A J2\{n\Fl\0)\' S{u; - (E^ - Eo)) , (44) 

n 

where the composite fermion operator is given as 



ImT^(t^) = 2sin\5/-K) 



(|44(l can be determined from the 
Tff\YQy-Q § denotes the 



The constant A in Eq 
condition 

phase shift associated with the fixed point spectrum, and 
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FIG. 9: (color online) The imaginary part of the fermions' T- 
matrix for jz ~ 0.15, j± — 0.15 and for different values of the 
bosonic coupling. For decoupled bosonic heat bath, the usual 
Kondo resonance shows up. As a is increased, the width of 
the Kondo resonance is decreasing and at the critical bosonic 
coupling (qc — 0.304) it disappears and the electrons become 
decoupled from the spin. 



is (5 = 7r/2 at the quantum fluctuation dominated Kondo 
fixed point. 

The results are shown in FiglHI These results show 
that in the fermionic phase the dot-lead junction opens 
up, and the transmission between the box and the lead 
becomes perfect at T = temperature. It is adequate to 
mention here that the T-matrix computed here is much 
less accurate than that of the one channel Kondo model. 
The reason is that due to the much bigger Hilbert space 
we lose a substantial part of the spectral weight. 



V. CONCLUSIONS 

In the present paper we studied the dissipation-induced 
phase transition at the degeneracy point of a single elec- 
tron box, discovered by Le Hur, using non-perturbative 
methods such as bosonization and numerical renormal- 
ization group (NRG). Throughout this paper we focused 
to the limit where the level spacing of the box is neg- 
ligible. Our calculations confirm the picture of Ref. |3 
For small dissipation the dissipative coupling is irrelevant 
and charging steps of the isolated dot become rounded. 
However, above a critical value the dissipative coupling 
caused by fluctuations of the electromagnetic field domi- 
nates, and suppresses quantum fluctuations of the charge 
on the dot. As a results, the charging steps are restored 
for large enough dissipation. 

Our non-perturbative analysis reveals that - as already 
found by the perturbative analysis of Ref. 9 - this phase 
transition is of Kosterlitz-Thouless type. Consequently, 
the Kondo energy, where quantum fluctuations start to 
dominate goes to zero exponentially fast as a approaches 
its critical value ac. 

This phase transition also shows up in the conductance 
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of the junction, which has a jump. This can be seen as 
foUows. For a > Uc the scahng equations imply that j± 
scales to 0. Since the linear conductance between the lead 
and the box is proportional to j\ ^ t^, this implies that 
the liner conductance for a > ac is identically zero. For 
a < Q!c, on the other hand, the system flows to the Kondo 
(two-channel Kondo) fixed point, where the conductance 
is 5 = 1 both in the spinless and in the spinful case. As a 
result, the T — Q temperature conductance has a jump: 



g{T = 0, a) = ^0{ac - a) 



(45) 



We can also determine the conductance at the critical 
value of a as a function of temperature from our non- 
perturbative equations. From the scaling equations we 
obtain that at a — ac the coupling j± scales to zero 
approximately as 



1 



ln(r*/T) 



(a = ac) 



(46) 



where T* is a dynamical scale. As a result, the box-lead 
conductance scales to zero as 



1 



In^ T* IT 



[a = etc) 



(47) 



The NRG method makes us possible to determine the 
phase boundary for the case of spinless Fermions between 
the two phases non-perturbatively. For small values of 



the tunneling we find from the scaling analysis that 



(48) 



which agrees very well with the NRG results. This con- 
dition translates to the following condition for the device 



4C2 Rq 



':V9 



(49) 



where g is the dimensionless conductance of the junction, 
Cg is the gate capacitance, C is the total capacitance of 
the island, and Rg is the resistivity of the gate electrode. 
In other words, the above phase transition can be reached 
by increasing the resistance of the gate electrode or by 
changing the box-lead conductance. 

Let us finally make a remark on the finite level spacing 
A. All considerations above hold only in the limit of van- 
ishing level spacing. In reality, however, this level spacing 
is finite, and plays a role similar to a finite temperature. 
As a consequence, quantum fluctuations remain always 
finite at the degeneracy point, and it is not clear if a real 
phase transition occurs even in the presence of a finite 
level spacing. 
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